Bragg Spectroscopy of Cold Atomic Fermi Gases 
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We propose a Bragg spectroscopy experiment to measure the onset of superfluid pairing in ultra- 
cold trapped Fermi gases. In particular, we study two component Fermi gases in the weak coupling 
BCS and BEC limits as well as in the strong coupling unitarity limit. The low temperature Bragg 
spectrum exhibits a gap directly related to the pair-breaking energy. Furthermore, the Bragg spec- 
trum has a large maximum just below the critical temperature when the gas is superfluid in the 
BCS limit. In the unitarity regime, we show how the pseudogap in the normal phase leads to a 
significant suppression of the low frequency Bragg spectrum. 
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I. INTRODUCTION 

Bragg spectroscopy, using two laser beams, has proven 
to be a very successful tool for probing the structure fac- 
tor of trapped atomic BEC's [lj. In this paper, we exam- 
ine the analogous experiment on two-component Fermi 
gases in both the normal and superfluid phases. We 
calculate the Bragg scattering rate for varying interac- 
tion strengths in both the weak coupling BCS and BEC 
limits and in the strong coupling unitarity limit. Previ- 
ously, we proposed measuring the inelastic (Stokes/anti- 
Stokes) scattering of an off-resonant laser beam on a two- 
component Fermi gas Q- In the BCS limit, the inten- 
sity of the scattering in this Stokes experiment exhibits a 
large maximum just below the critical temperature, T c , 
at which the gas becomes superfluid. This effect is the 
light scattering analog of the Hebel-Slichter effect in con- 
ventional superconductors, a hallmark experimental test 
of BCS theory y|. Here, we demonstrate that the same 
effect can be observed in a Bragg scattering experiment. 
We furthermore show that the pairing induced by the 
interactions can be readily detected in the Bragg spec- 
trum. In the strong coupling regime, we discuss how the 
presence of a pseudogap in the normal phase has signifi- 
cant effects on the Bragg spectrum. The advantage of the 
Bragg experiment as compared to the Stokes experiment 
is that one varies in a controlled way the energy and mo- 
mentum imparted to the atoms. Compared with recent 
radio frequency (rf) experiments which involve a third 
unoccupied hyperfine state £| , the Bragg experiment in- 
volves only the two hyperfine states already present in 
the trap. The Bragg experiment, avoiding complications 
due to non-trivial interaction effects between the third 
hyperfine state present in the rf experiment [5| and the 
trapped gas, is thus simpler to interpret. 
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FIG. 1: Schematic atomic level diagram for stimulated 
scattering leading to a "spin" flip. 



&b = 1 in this paper) above the state | |), and that the 
population of the two states is equal, = N± = N/2. 
The two atomic states interact with a coupling which can 
be characterized by a scattering length a in a vacuum. 
We ignore the interaction between atoms in the same 
hyperfine state since s-wave scattering is suppressed be- 
tween two identical fermions. 

Bragg scattering is realized by a pair of laser beams 
with wave vectors ki and 1<2 and frequencies u>i and ioi 
illuminating the gas. An atom absorbing a photon from 
beam 1 and emitting a photon in beam 2 changes its 
energy-momentum by w\ — tt»2,ki — \s.2- Bragg scatter- 
ing can be used to probe both the density and spin cor- 
relation functions. Pairing is predicted to have several 
effects on the frequency and momentum dependence of 
these correlation functions for T = 

IS Ell- To probe 

the spin-flip correlation function, one tunes the two laser 
frequencies so that lo\ — 0J2 ~ ojm as illustrated in Fig.^ 
If the detuning, 5, from the optical transition to or 
from the upper level \e) is large, population of the up- 
per level can be neglected, and the effective light matter 
Hamiltonian describing spin-flip Bragg scattering can be 
written as 



II. BRAGG SCATTERING 

We consider a homogeneous gas of fermionic atoms oc- 
cupying two hyperfine states, denoted by | f) and | [). 
We assume that the state | [) has energy (% = 1 and 



rf 3 , 



h.c. 



(1) 



where X(r, t) 



-^ 2)t l(r), and Z(r) describes the 
spatial profile of the laser beams and the dependence 
on the atomic states involved and their coupling to the 
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electromagnetic field We neglect finite waist effects 
and take Z(r) = Xe lq r with q = — k 2 . 

To measure the effects of pairing we focus on the spin- 
flip rate (N\) induced by the Bragg lasers, where N± is 
here the number of particles in the level | J,), and (. . .) 
denotes a thermal average. The Heisenberg equation of 
motion for the operator N± is 



[N l7 H E 



d A r 



h.c. 



(2) 

For sufficiently weak laser beams, (N\) can by calculated 
by linear response theory: 




-2Z 2 ImX>(q,w + w/), 



(3) 



with uj — io\ — LU2 — Whf. Bragg scattering thus probes 
the Fourier transform of the spin-flip correlation func- 
tion ([^(r,t)ipi(r,t),ipl(r',t')ipt(r', t')]). The response 
of the gas is usually detected experimentally by measur- 
ing the total momentum imparted on the gas |l| ; the rate 
of momentum transfer is straightforwardly related to the 
scattering rate by 

P = q^l (4) 

where P = P j + Pj is the total momentum of the gas in 
the two hyperfine states. 



III. WEAK COUPLING BCS LIMIT 

In the weak coupling limit a — > 0~, the gas can be de- 
scribed by BCS theory. The correlation function in J3J is 
then a sum of contributions from normal and anomalous 
Green's functions ^(j- We obtain after some algebra, 

V o(q,z) = [ 7^ [(««' + vv' f 



1 



[uv — vu ) 



(2tt) 3 
(1-/ 



/') 



E-E< 



f-f 
z + E-E 

(f + r-iy 

z + E + E' 



(5) 



with z = lu + irj. Here, u = u k and u' = itk+qj etc., with 
u\ = (l + £k/E k )/2 and v\ = 1 — u\ as usual in BCS the- 
ory. The quasiparticle energies are E — Ek = + A 2 
and E' = -Ek+q with = ek — /i, where e k = k 2 /2m, 
fi is the chemical potential, and A the BCS pairing gap. 
The Fermi functions are / = ,f(E) and /' = f(E'), with 



,x/T 



l) . The first term of (J5J, describes 



scattering of quasiparticles while the last two terms de- 
scribe the creation and annihilation of two quasiparticles 
- pair breaking processes. 

For an ideal gas in the normal phase, 10) yields in the 
limit u> -C ep, T <C Tp (Tp is the Fermi temperature), 
and u> <C qkp/m, 

I m V ( q ,u) = ---^, (6) 
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FIG. 2: The scattering rate —lmD(q,Lj) in units of the scat- 
tering rate in the normal phase vs. T for representative u>. 
(color online). 



where ep is the Fermi energy, n — kp/3ir 2 is the total 
density of the gas, and kp is the Fermi momentum. 

The spin flip induced by Bragg scattering can excite 
collective spin waves. These modes are undamped in the 
normal phase for an attractive interaction [Tlj . For a re- 
pulsive interaction in the normal phase, the spin waves 
are damped and they change the spectral shape of the 
particle-hole continuum. The single particle approxima- 
tion JSJ does not include such collective modes. In order 
to include them, we sum bubbles to obtain the random 
phase (RPA) approximation [To| . 



2?rpa 



1 + TV ' 



with r = 4ira/m the scattering matrix for fcpa <C 1. 
Note the + sign in the denominator due to the spin flip. 



A. Results 

We first consider the temperature dependence of the 
scattering rate. An attractive feature of the Stokes ex- 
periment is that one can observe a strikingly large max- 
imum in the scattered signal just below the critical tem- 
perature T c 0- We now demonstrate that the same 
effect is present for the Bragg scattering experiment in 
the experimentally realizable parameter range, u> <C A, 
U) < qkp/m, and q < 2kp. 

Figure |21 shows the scattering rate J3Jl as a function of 
T in the weak coupling limit, in units of the normal phase 
rate, which is constant for T <^ Tp. For illustration, we 
have taken fcpa = —0.3 and q = kp. Solving the BCS 
equations as a function of T, we obtain the scattering 
rate from Eqs. JSJ and Q. The critical temperature is 
T c = 0.0033eF- We plot the results for frequency differ- 
ences between the two laser beams lo = 0.031T C , 0.92T C , 
and 6.1T C . The numerical normal phase results agree 
with 10 , confirming the numerical accuracy of the calcu- 
lations. 
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As Fig. [5] shows, the scattering rate in the super- 
fluid state has a large maximum below T c for u < T c . 
This maximum is due to a large quasiparticle scatter- 
ing rate arising from the increased density of states at 
the Fermi energy in the superfluid phase just below T c , 
N(E) = (p F /2n 2 )E k /£ k . For u < T c , both scattering 
states k and k + q in Eq. JSJ can be in the region around 
the Fermi surface of a large density of states yielding an 
increased scattering rate. The maximum in the scatter- 
ing rate depends crucially on the coherence factors for 
quasiparticle scattering in (jSJ) adding, viz., uu' + vv' ~ 1, 
for a spin flip. For Bragg scattering into the same hy- 
pcrfine state, the coherence factor is uu' — vv' ~ at 
the Fermi surface and the increased density of states 
would not lead to a maximum in the scattering rate. For 
T — > 0, the scattering rate is suppressed by the increasing 
gap A(T) since the density of quasiparticles available for 
scattering scales as exp(— A/T). For u> < 2A the lasers 
cannot break Cooper pairs. 

The physical reason for the coherence factors adding 
(uu' + vv' ~ 1) for spin-flip scattering is that the quasi- 
particles carry the same spin in the superfluid as in the 
normal phase. On the other hand, in scattering without 
spin flip, which couples to the particle density, the parti- 
cle number carried by the quasiparticles is suppressed at 
the Fermi surface, since a quasiparticle is a superposition 
of a hole and a particle. In this case the coherence factors 
subtract (to' — vv'). 

For increasing to, the peak in the rate below T c de- 
creases and eventually disappears for u) ~ T c . From 
Fig. |2 we see that the scattering rate in the superfluid 
state is always smaller than in the normal state when 
w ^ T c (green dashed and black dash-dotted lines). This 
is because both scattering states cannot be in the region 
around the Fermi surface where the single particle density 
of states is increased. Bragg scattering does not probe 
the increased density of states for the superfluid state for 
higher frequency as effectively as for co <C T c , and the 
resulting rate is decreased compared to the normal state. 
Note that the kink in scattering rate just below T c for 
u> = 0.92T C , at which point 2A(T) = to, is physical. Be- 
low this temperature, the gap is too large to allow the 
creation of two quasiparticles and the second term in (0 
vanishes. 

Finally, Fig. [21 for cl> = 6.1T C illustrates the case of 
large energy transfer, u> > 2A(T = 0). Contrary to the 
case oj < 2A(T = 0), the rate is finite for T — as the 
energy transfer is large enought to break pairs and the 
second term in (JSJ is finite. For oj/T c — > oo, the scattering 
rate approaches that of the normal state, as expected. 

We conclude that there is a peak in the scattering rate 
for T <T C and uu <C T c . Also, the peak disappears unless 
lo < qkp/m and q < 2k p; otherwise, both scattering 
states cannot be in the region of increased density of 
states. 

We now consider the frequency dependence of the 
Bragg scattering rate. Figure [3] shows the rate from Q 
as a function of u> for the gas in the weak coupling limit 
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FIG. 3: The scattering rate —lmD(q, u>) in units of n/2eF as a 
function of u) for fc^o = —0.3 and various temperatures (color 
online) . 



with kpa = —0.3 and q/kp = 0.0035. The scattering 
rate is depicted for various temperatures below and just 
above T c where the gas is normal. For comparison, we 
also plot the result of the single particle approximation 
in the normal phase. 

The large peak in the normal phase response just above 
the single particle continuum arises from the spin wave 
present above the particle-hole continuum for an attrac- 
tive interaction. The finite width of the peak is an arti- 
fact of our having included a small imaginary part in the 
frequency u> for computational reasons. We see from Fig. 
|3 that when the gas becomes superfluid, the spin wave 
vanishes. This is because only potential flow, i.e., den- 
sity fluctuation modes, is possible in the superfluid [l2T |. 
Other modes such as spin waves require a normal compo- 
nent to be present. In a trapped system this effect leads 
to the damping of the spin-dipole mode in the superfluid 
phase |13j- In the superfluid phase, the particle-hole spec- 
trum disappears and is replaced by a low frequency gap, 
which reflects the pair breaking energy described by the 
second term in (0). At zero temperature the gas does 
not respont for ui < 2 A, the minimum energy required to 
break a pair. 

The pairing gap can thus be detected directly as a 
gap of 2A in the Bragg spectrum. By contrast, in the 
rf experiment of the shift in the spectrum scales as 
A 2 /cf, since one transfers particles to a third hypcrfinc 
state which in general interacts differently with the two 
paired hyperfine states (§(]• 



IV. WEAK COUPLING BEC LIMIT 

We now consider Bragg scattering in the weak cou- 
pling BEC regime kpa — > + . In this limit, the molecules 
are tightly bound with a radius much smaller than the 
average interparticle distance, and are weakly interact- 
ing both with each other and with unbound atoms |l4j . 
The system can therefore to a good approximation be 
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regarded as an ideal gas of point molecules and atoms in 
thermal equilibrium. The molecules Bose condense be- 
low the critical temperature T c = 0.2187V. The Bragg 
scattering induced by Q then either breaks a molecule 
by flipping the spin of a bound | atom, thereby creating 
two free j atoms, or it simply flips the spin of a free | 
atom. The Fermi Golden Rule result for Bragg scatter- 
ing rate into the free atomic J. states for such a system is 
then a sum of three terms, 



hi=i (r at 



r 



BEC 



The scattering from the unbound atoms is given by 
d 3 k 



(27T>- 



r Cf-/')27r5(w 



(8) 



(9) 



where e = e k , e> = e k+q , / = /(£), and /' = /(£'). The 
rate from a BEC of molecules with density n is 



Tbec — no 



d 3 k 



\M\ 2 (l~f-f')2ir5(Lu + E b - 



(10) 

with E b = —h 2 /ma 2 the energy of a molecule with zero 
momentum relative to the bottom of the free particle con- 
tinuum of a pair of f and J, atoms, and M = M(k, q) the 
matrix element of Q (excluding X) between a molecule 
at rest and a pair of free J, atoms with momenta k + q 
and — k. The scattering rate from thermal molecules is 



Mol 



d 3 p f d 3 k 



(2tt) 3 J (2tt) 3 



[n B (e mo l)(l - / - /') - ff]2*6(w + E b + ^--e~ e'), 

4m 

(11) 

with ns(x) — [e x l T — 1) _1 the Bose distribution, and 
e m oi(p) = E b + p 2 /4m — 2[i the energy of molecules with 
momentum p. The molecule breaks into two j atoms 
with momenta k + q + p/2 and p/2 — k and energies 
e = e p /2-k an d e ' = e k+q+p/2- The matrix element in 
(flip is, from Galilean invariance, independent of p. 

To proceed, we model the molecule-free atom pair ma- 
trix element M — M(k, q) as follows. We write state of 
a molecule at rest as \cj>) = X)k ^fc c kT c -kjjO) > w ^h |0) the 
vacuum and c fc(7 creating an atom with momentum k in 
hyperfine state a. The Fourier transform of the molecu- 
lar wave function is <f>k . The spin- flip process illustrated 
in Fig. ^ breaking the molecule and adding momentum 
q creating two free J. atoms with momenta k + q and 



-k corresponds to f^+qiCktl^) — ^fe c 



t 

k+q| 



JO). How- 



ever, the process c_ k ^c_k 



-qTl 



c -kii 



|k+q| L k+ql L -kll 



connects to the same final state, and the two matrix el- 
ement must be added coherently in the Golden Rule ex- 
pression. We thus find the matrix element |M(k, q)| 2 = 
\4>k ~ <^jk+q|| 2 /2 (where the factor 1/2 avoids double 
counting) . 
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FIG. 4: Contributions to the total atom density (in units of 
n) from the molecular BEC, the thermal molecules, and the 
free atoms. 



The BCS wave function correctly describes a conden- 
sate of non-overlapping bosons (Cooper pairs) with bind- 
ing energy E b = —h 2 /ma 2 in the BEC limit a — * + . 
To connect this wave function to the above analysis, we 
note that in the BCS expression the Bragg scattering 
(J5J, it is the second term that describes the breaking 
of Cooper pairs. Thus, BCS theory yields |M(k, q)| 2 = 
(wfeUk+q - ffeUk+q) 2 /2. In the BEC limit, we have Uk — > 1 
and v k -> A/(2e k + \E b \) with A -> '{my/a). BCS 

theory therefore yields <pk for the Fourier transform 

of the molecular wave function and 



Vke oc 



,-r/o 



lar 



(12) 



in the BEC limit. This is the asymptotic bound state 
wave function for a potential with scattering length a as 
expected. We therefore use the BEC limit of the BCS 
result for the matrix element |M(k, q)| 2 . Note that we 
neglect any closed channel components deep in the BEC 
regime. 



A. Results 

The populations of the molecular and atomic states are 
found by solving the atom-molecule equilibrium problem; 
the density of the thermal molecules, the molecular BEC 
and the thermal atoms are shown in Fig.0]for varying T. 

In Fig.0 we plot the scattering rate as a function of id 
obtained from JSJ for a molecule energy E b = — 22.22ei? 
(corresponding to kpa = 0.3) and for momentum trans- 
fer q = kp. The scattering rate is shown in Fig. for 
a pure molecular BEC at T = (a) and for thermal 
molecules at T — 0.5Tp above the BEC critical tem- 
perature T c = 0.218T F . For T = 0, the rate is given 
by p0(l and the threshold for the scattering is given 
by lu > —E b + q 2 /Am - the minimum energy required 
to break a p = molecule and create two J, atoms 
with momenta q/2 and — q/2. The rate from the ther- 
mal molecules at T = 0.5Tp obtained from I jlljl does 
not differ significantly from that of the molecular BEC 
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FIG. 5: Bragg scattering rate in units of nl 2 n/tF: from the 
molecular BEC, the thermal molecules, and the free atoms. 



at T = 0. This is because both temperatures are so 
low that there are essentially no atoms present, as one 
sees in Fig. QJ The only effect of finite temperature is 
then a small shift u> > min(— E% + q 2 /Am — qp/2m) ~ 
— Eb + q 2 /Am — q^fksT /m in the threshold for scatter- 
ing, from (|ll|l . Thus, the onset of Bose-Einstein conden- 
sation of molecules has only minor effects on the Bragg 
scattering rate since \Eb\ ^> ksTc in the BEC limit. 

In Fig. 03(b) we plot the scattering rate for T = 2.5Tp, 
at which most of the atoms are still bound as molecules 
but some free atoms are present, as can be seen from Fig. 
0J The Bragg scattering consists of a signal down to u> = 
from the free atoms, given by Q, and a signal from the 
thermal molecules for ui > — Ef,, given by Qllfl. Note that 
even though the atomic signal at small to is much larger 
than the molecular signal at high u>, the contribution to 
the f-sum rule from the molecules is much larger than 
that from the free atoms. For higher temperatures, the 
molecular signal of course disappears with the vanishing 
molecule density, and the Bragg scattering is given by 
that of the free atoms only. 



V. STRONG COUPLING UNITARITY LIMIT 

The scattering length diverges, kp\a\ — > oo, at the Fes- 
hbach resonance. In this unitarity limit, the interactions 
are strong and considerably modify the properties of the 
gas. An important characteristic of the Fermi gas at uni- 
tarity is the presence of a pseudogap in the normal phase 
close to T c ; i.e., the suppression of the single particle DOS 
around the Fermi surface 0] . A pseudogap has been ob- 
served experimentally for high T c superconductors where 
it is presently subject to intense investigations [16j. Since 
a gap in the spectrum of single particle excitations nat- 
urally leads to a gap in the low temperature Bragg spec- 
trum (Figs.EEJ, we now examine whether the pseudogap 
can be detected in Bragg spectroscopy. 

We calculate the Bragg spectrum at unitarity using a 
minimal many body theory that includes the correct two- 
particle physics leading to the resonance. For the single 
channel problem, such a theory including pair fluctua- 



tions in the ladder approximation was originally devel- 
oped many years ago fl7| and later examined in great 
detail by several authors |14lT5| . This approach was then 
extended to the multichannel case relevant for atoms in- 
teracting via a Feshbach resonance , tr eating the bound 
state as a point boson [l8L Il9l I20L l2l| , and taking into 
account that the bound state is a composite two-fermion 
object IH. 

The atoms are described by the atomic propagator 
G(q, z)- 1 = G (q, z)" 1 -^, z) with G (q, z)- 1 = 
The atom self-energy is given in the ladder approxima- 
tion by 

E(<?,w n ) = /r 1 Tr[r(K,p,p,c Jn + Wm )Go(A : ,w m )] (13) 

where K = q + k, p = (q — k)/2, the trace denotes a 
sum over Matsubara frequencies u> m = i(2m + 1)ttT and 
integration over k, and T is the many-body scattering 
matrix depending on the center of mass K and relative 
momenta p of the two scattering particles [lj}. -"- n the 
limit of a broad resonance, the effective interaction me- 
diated by the Feshbach molecule can be regarded as in- 
stantaneous with a resonant scattering length \a\ — > oo 
at unitarity plj . Neglecting the frequency and momen- 
tum dependence of the effective interaction coming from 
the molecular state, we recover the single channel theory 
scattering matrix 



T(K,z)- 1 = — -U(K, 



Ana 



The pair propagator is, as usual 
U(K,z)-- 



d 3 k (1-f-f m 



(14) 



(15) 



Here, £ = £ K /2+k, ? = £ K /2-k, / = f(0, and /' = /(£')• 
The structure of Eqs. l(T5 )l -ljT5 )l is indicated in Fig. 0(a). 

To obtain the Bragg spectrum, we need to evaluate the 
atomic propagator at frequency to + irj. (From here on 
we do not write the infinitesimal positive imaginary part 
of the frequency explicitly.) From (| 1 3|) . we obtain 

X(q,w) = J ^ j ^-[f(e)A (k,e)T(k\e'y 

n B (e')A m (k',e')G a (k,e)}. (16) 

Here k' = k + q, e' = e + ui, and Ao(k,u>) — 
— 2ImGo(fc, u>) and A m (k, w) = — 2Imr(fc, u>) are the spec- 
tral functions for non-interacting atoms and the pair 
propagator, respectively. We evaluate the real and imag- 
inary parts of E directly from (|16f) . We also calculate 
the real part from the imaginary part using a Kramers- 
Kronig relation, and find that the two calculations agree 
to a high accuracy, confirming the consistency of the nu- 
merical calculations. 

In order to evaluate medium effects, we need to deter- 
mine T c and the chemical potential, fi. We determine fj, 
by calculating the thermodynamic potential Q at unitar- 
ity in the ladder approximation. The correction Ail to 
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FIG. 6: (a) The atomic propagator in the ladder approxima- 
tion for a broad resonance where the Feshbach interaction is 
effectively instantaneous at unitarity; (b) the spin-flip corre- 
lation function in a conserving ladder approximation. 




FIG. 7: The atom spectral function A(k, uf) at unitarity and 
T = T c for various momenta. 



the ideal gas result can be expressed as an integral over 
the phase of the on-shell scattering matrix T [23| : 



d*Kd 3 k mm 

(2tt) 6 ImII(K,£ + £' 



■9(K,C + e) (17) 



with £, £' defined below l|15|l . The angle 9(ui) is defined 
by writing T(lj) = \T\e l6 . Equation (|1 T|) reduces to the 
virial expansion for Q at high temperatures. The chemi- 
cal potential \x is then determined by —dVl/d^ = n. 

In order to calculate the spin-flip correlation function 
consistently, one needs in principle to choose a func- 
tional $[G], generate the self-energy as E = 5&/6G, and 
then generate the correlation functions from the kernel 
<5E/(5G |24|. The diagrams relevant for the ladder ap- 
proximation are given in Fig. [B| (b) where all propagators 
are full Green's functions, G, including those in T. The 
conserving spin-flip correlation function thus obeys an 
integral equation with T and the full G's. The imaginary 
part of the first term in this equation is 



T ^ , \ 1 f d 3 k f°° de „ , .. 



-A (18) 



with A = A(k,uj), A' = A(k',u>), and k' and £' defined 
as in (|16[) . The spectral function A(k,cv) = — 2ImG(fc, w) 
describes the excitations of the interacting atoms. 



A. Results 

We calculate the Bragg scattering rate at the critical 
temperature T c where we expect the pseudogap effects to 
be most pronounced. From the thermodynamic potential 
(|17|l . we obtain /i(T) and can calculate the critical tem- 
perature from the Thouless criterion r _1 (0, 0, T c , /i) = 0. 
We then evaluate the atom self-energy from i|16[l and from 
this obtain the spectral function. 
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FIG. 8: (a) The positions E b (k) and E p (k) of the hole and 
particle peaks of A{k,u>) as a function of k at unitarity and 
for T =T C . The solid line is the free particle energy k 2 /2m — 
fi. (b) The corresponding spectral weights of the particle 
and hole peaks. The dashed lines are fits to simple analytic 
functions. 



To understand the Bragg scattering rate at unitarity, 
we first analyze the effects of the interactions on the sin- 
gle particle spectrum, described by the spectral function, 
A(k,oj). In Fig. [7| we plot A(k,u>) at T c for various 
momenta, measured in units of fc M = ^/2m/i. We have 
chosen parameters corresponding to a resonant interac- 
tion with \kpa\ — 11.8 ^> 1 and a negligible effective 
range. For this set of parameters, we find T c sa 0.26Tp 
with fj,(T c ) m 0.45ei? in good agreement with other BEC- 
BCS crossover results based on a similar approxima- 
tions |l8|, |25j ■ The spectral functions obey the sum rule 
J A(uj)duj /27r = 1 to a very good approximation, pro- 
viding another important check of the numerical calcula- 
tions. 

In Fig. [7| we clearly see the double peak structure of 
A(k 7 to) with the spectral weight at u> = significantly 
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FIG. 9: — lmDo(q, to) at unitarity and T — T c in units of 
n/2eF as a function of for(color online). 

suppressed for all k. This suppression is the characteristic 
feature of the pseudogap. We refer to the peak at u < 
as the hole peak at frequency —E^, and the peak at lo > 
as the particle peak at E p . The positions of the two peaks 
found by fitting to a double Lorentzian (normalized such 
that its total spectral weight is 2ir) are plotted in Fig. |S| 
(a). 

We see that the distance of the maxima to lo = 
reaches a minimum: min£ ; h(fc) ~ mmE p (fc) ~ 0.6ef at 
k r~j 1.2k fj,, providing a qualitative value for the pseu- 
dogap. The spectral weight of the hole peak decreases 
and that of the particle peak increases with increasing k. 
This behavior is illustrated in Fig. [S] (b) where we plot 
the spectral weights U 2 of the particle peak and the hole 
peak 1 — U 2 obtained by fitting to the double Lorentzian. 
For large k/k^, the spectral function approaches that of 
an ideal gas with a single peak at lo — k 2 /2m — ji as 
expected. The behavior of the spectral function at uni- 
tarity as described here is in general agreement with the 
analysis of Ref . . The spectral function behaves qual- 
itatively as the BCS spectral function 

A BCS {k,Lo) = v 2 2it5(lo + E k ) + u 2 2ir5{uj ~ E k ), (19) 

with two peaks located around to — 0. There are impor- 
tant differences from BCS theory though. The two peaks 
for the spectral function in the unitarity limit at T c have 
large widths (which are unequal for the particle and hole 
peaks) and the spectral density is suppressed but finite 
at to = 0. The pseudogap is a suppression of spectral 
density and not a real gap. Also, the positions of the two 
peaks are not symmetric around to = and the minimum 
gap is not at k = k^ but slightly above, as can be seen 
from Fig. [S] 

We now ask whether the pseudogap behavior of the 
spectral function at unitarity has consequences for the 
Bragg spectrum. In Fig. [51 we plot lmT>o(q,Lo) for 
q — k^/2 and T — T c . The x's are obtained from JTSJ 
using the spectral function calculated from ■ We eval- 



uate the (computationally demanding) multidimensional 
integrations using a Monte Carlo routine. The solid line 
in Fig. [5] is obtained by using the double Lorentzian fit 
for the spectral function in (|18|) . Here, the po- 
sition of the hole Eh(k) and particle E p (k) peaks, the 
width and the weight U 2 (k) are calculated as functions 
of momentum k using approximate functions obtained 
from fitting A(k,u>) to the double Lorentzian for var- 
ious k. The approximate functions for Eh{k), E p (k), 
and U 2 (k) are given as dashed lines in Fig. [S] Once 
the approximate double Lorentzian form A^(k,Lo) is ob- 
tained, the integration in (|18H is much faster than when 
using the spectral function A(k, lo) obtained from the full 
ladder calculation. We see from Fig. that the double 
Lorentzian model A^k^uj) yields results for ImX>o(<?, a;) 
in good agreement with the full numerical calculation 
based on A(k,u). From this, we conclude that the pseu- 
dogap spectrum close to T c can be well described by a 
simple double peak spectral function. 

Comparing to the ideal gas scattering rate for the 
same density and temperature, also plotted in Fig. El 
we see that interactions suppress the spectral weight of 
lmT>o(q,w) at low u> and push it to higher frequencies. 
This is a direct effect of the pseudogap. As opposed to 
the real gap in the T — Bragg spectrum, the spectral 
density is however finite at small lo. This is as expected 
since the single particle spectral density is suppressed but 
finite for lo = (Fig. EJ). Also, T c ~ 0.26e^ is compara- 
ble to the pseudogap ~ 0.6ei? so thermal excitations are 
present. From l|18f) . we see that that both effects make 
ImT>Q(q,Lo) finite for small lo. 

The approximation T>o(q, lo) to the spin-flip correlation 
function is not conserving and does therefore not obey 
the f-sum rule. To obtain a conserving approximation, we 
need to solve the integral equation in Fig.El(b). This will 
change the quantitative form for T>(q, lo). The solution is 
however an iteration with the Y matrix of T>o(q,u>). We 
therefore expect Iml?(q, lo) to have the same qualitative 
characteristics as ImT> (q, lo), i.e., a suppressed but finite 
spectral weight at low lo. Solving the fully self-consistent 
equation for the spin-flip correlation function is beyond 
the scope of this paper [2(j . 

In the BCS limit, we see that the emergence of a gap 
in the single particle spectrum leads to a maximum in 
the Bragg scattering rate just below T c (Fig. |2J. One 
could therefore expect a similar maximum to appear at 
unitarity at a temperature T* defined as the temperature 
where the pseudogap emerges. However, the presence of 
a maximum requires the spectral function to have well- 
defined sharp quasiparticle peaks like in the BCS case 
l[T5)l. The pseudogap spectral function at unitarity (Fig. 
EJ on the other hand has peaks of width ~ C(e_F )■ There 
will therefore be no maximum in the Bragg scattering 
rate for T ~ T* due to the emergence of the pseudogap 

We conclude that the pseudogap present in the normal 
phase close to T c at unitarity leads to a significant sup- 
pression of the low frequency Bragg scattering rate. The 
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rate remains however finite since the pseudogap is only 
a suppression of the single particle spectral density, and 
because T c is comparable to the value of the pseudogap. 
The emergence of a pseudogap with lowering T does not 
give rise to a maximum in the scattering rate, as opposed 
to the BCS case. This is because the quasiparticles are 
strongly damped at unitarity. 

VI. EXPERIMENTS 

We briefly comment on the possible experimental ob- 
servation of the effects described in this paper. A typical 
Fermi energy for trapped atomic gases is <~ O(10 4 )Hz. 
Thus, we require a frequency resolution of the Bragg ex- 
periment <C 10 4 Hz. Likewise, a typical Fermi momentum 
of the gas is &f ~ O(10 6 )m _1 which is the same order of 
magnitude as the wave numbers for optical lasers. This 
means that the resulting angle between two laser beams 



needed to obtain a momentum transfer of ~ kp is rea- 
sonable. The results presented in this paper are for a 
homogenous system. From the Thomas-Fermi approxi- 
mation, we expect the presence of the trapping potential 
to smooth out features like the gap in the Bragg spec- 
trum. To avoid such smoothening, one could focus the 
laser beams on a smaller section of the cloud where the 
gas can be regarded as homogeneous. In this case, the 
results presented in the present paper become directly 
applicable. 
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